Self-motile colloidal particles: from directed propulsion to random walk 
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The motion of an artificial micro-scale swimmer that uses a chemical reaction catalyzed on its 
own surface to achieve autonomous propulsion is fully characterized experimentally. It is shown 
that at short times, it has a substantial component of directed motion, with a velocity that depends 
on the concentration of fuel molecules. At longer times, the motion reverts to a random walk with 
a substantially enhanced diffusion coefficient. Our results suggest strategies for designing artificial 
chemotactic systems. 
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The directed propulsion of small scale objects in wa- 
ter is problematic because of the combination of low 
Reynolds number and Brownian motion on these length 
scales lj . In order to achieve an artificial micro- or nano- 
scale swimmer that is able to propel itself in a purpose- 
ful way, one needs both a swimming strategy that works 
in the environment of low Reynolds number [2|, and a 
strategy for steering and directing the motion that can 
overcome the ubiquity of Brownian motion. Common 
bacteria, such as E. Coli, achieve propulsion by non-time- 
revcrsible motion of long flagella, and employ a "run and 
tumble" strategy to be able to swim towards or away 
from environmental stimuli 0]. 

One possibility for designing propulsion is to devise 
non-reciprocal deformation strategies that are simple 
enough to be realizable Recently, an interesting ex- 
ample of such robotic micro-swimmers has been made us- 
ing magnetic colloids attached by DNA-linkers, and con- 
trolled by an external oscillatory magnetic field An- 
other possibility is to take advantage of phoretic effects, 
where gradients of fields such as concentration, tempera- 
ture, electric field, etc, couple to the surface properties of 
particles to create slip velocity patterns that could lead 
to net propulsion [f|. 

A particularly appealing strategy for propelling small 
devices is to take advantage of chemical reactions Q. 
In a pioneering experiment, Paxton et al. observed au- 
tonomous motion of platinum/gold nanorods Q, and 
similar experiments have been performed by Fourneir- 
Bodoz et al who used gold/nickel nanorods [9( and Mano 
and Heller who used enzymes [10j. One simple strat- 
egy for converting chemical energy to mechanical work 
in such devices has been proposed by one of us and col- 
laborators [HI, [12]. I n this scheme, a spherical particle is 
considered with an asymmetric distribution of catalyst on 
its surface. If the chemical reaction so catalyzed produces 
more products than it has reactants, then the asymmetric 
distribution of reaction products propels the particle by a 



process of self-diffusiophoresis. The experimentally real- 
ized swimmers mentioned above are very similar in setup 
to the one proposed in Ref. 111 ]. However, it seems that 



the driving mechanism for the propulsion in these exper- 
iment is different as it depends strongly on the presence 
of a (bi- metal or bio-) electro-catalytic structure [l3| ■ It 
is suggested that the pair of electrochemical reactions at 
the two poles create a lateral electric field near the par- 
ticle surface (due to electron transfer) that could move 
the solvent via electro-osmosis 



Here we have realized the scheme proposed in Ref. [11 1 
experimentally by taking polystyrene spheres with nar- 
row size distributions with a diameter of 1.62//m, and 
coating one side of the spheres with platinum keeping 
the second half as the non-conducting polystyrene. The 
platinum catalyzes the reduction of a "fuel" of hydro- 
gen peroxide to oxygen and water, which produces more 
molecules of reaction product than consumed fuel. We 
fully characterize the motion of the artificial micro-scale 
swimmer using particle tracking, and probe the prop- 
erties of the motion as a function of hydrogen peroxide 
concentration. We show that at short times, the particles 
move predominantly in a directed way, with a velocity 
that depends on the concentration of the fuel molecules 
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FIG. 1: (color online). Trajectories over 25 sec for x5 par- 
ticles of the control (blank) and Platinum-coated particles in 
water, and varying solutions of hydrogen peroxide. 
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FIG. 2: (color). Mean squared displacement as a function of 
time interval for the Pt-coated sphere trajectories shown in 
Fig. [1] corresponding to different hydrogen peroxide concen- 
trations. 
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FIG. 3: (color online). Mean squared displacement as a func- 
tion of time for the inset trajectory (Pt-coated sphere in 10% 
H2O2) fitted to Eq. {TJ, which resembles a parabola (red line) 
at At <C tr and a straight line (green line) at At ^> tr. The 
fits yield v = 3.1 /ims -1 , D — 0.31 /im 2 s _1 , and tr = 3.9 s 
for this particular trajectory. 



with a Michaelis-Menten behavior. At longer times, the 
motion reverts to a random walk, in which runs of di- 
rected motion are interrupted by random changes of di- 
rection. We also extract the positional and the rotational 
diffusion coefficients and show that they are consistent 
with theoretical predictions, with the rotational diffu- 
sion featuring a moderate concentration dependence that 
could be attributed to directed rotational components in 
the motion. 

Polystyrene microspheres were produced as described 
by Fujii et al [lj]. The weight-average diameter of the 
microspheres was determined by dynamic light scatter- 
ing to be 1.62 ± 0.13/im, with a polydispersity index of 
1.008. To half-coat the spheres Hf in platinum, a di- 
lute suspension (0.05 wt%) of the spheres in isopropanol 
(HPLC grade: Fischer) was drawn over cleaned glass 
microscope slides to produce a dilute monolayer of the 
spheres (~ 10 5 spheres-cm -2 ). A thin layer of Platinum 
(99.99% Agar Scientific) was evaporated onto the micro- 
sphere coated slides resulting in a layer of Platinum 5.5 
nm thick on one side of the spheres. The half-platinum 
coated spheres were detached from the glass slide using 
a sheet of PTFE (200 /im thick - Goodfellows) acting 
as a blade and the detached spheres resuspended in 2 
ml of distilled water (18.2 MQ cm -1 - Purite). Hydro- 
gen peroxide solutions were prepared by subsequent di- 
lutions of 30 % Hydrogen Peroxide solution (Perdrogen - 
Riedel-de Haen) . For a given solution of hydrogen perox- 
ide, 100 /il of the microsphere containing stock solution 
was mixed with 2 ml of the appropriate solution. A cu- 
vette (1 mm optical path length - Hcllma, cleaned with 
Piranha Etch) was rinsed three times with the solution 
before filling again and sealed using a PTFE stopper. 
Particle tracking was achieved using a Nikon ME600 op- 
tical microscope mounted on an isolation table and fitted 
with a Pixelink PL-A742 machine vision camera, a x20 
objective (Nikon), and inverted illumination. A movie 



with a field of view of approximately 120 /im x 120 /im 
of 3000 frames at a rate of 38 frames per second was 
recorded for each particle. For each concentration of hy- 
drogen peroxide approximately 20 separate particles were 
tracked and analyzed using a Labview (National Instru- 
ments) script providing particle trajectory and timing 
(time, x/pm, and y//im). 

Figure [1] shows particle traces for a Pt-coated bead as 
well as a control polystyrene bead. While the control 
particle undergoes a characteristic Brownian motion, the 
Pt-coated particle shows a systematic enhancement of 
a directional component in its motion, as the concen- 
tration of hydrogen peroxide increases [16j | . The motion 
can be analyzed more quantitatively using particle track- 
ing. From the records of the particle trajectory, we cal- 
culate the average value of the squared displacement as 
a function of time. For a purely Brownian particle of ra- 
dius R, the squared displacement is linear in time with 
the slope controlled by the particle diffusion coefficient 
D = kBT/(6itr]R), where k^T is the thermal energy, 
r\ is the viscosity of water. The particle will also un- 
dergo rotational diffusion with a characteristic (inverse) 
time scale TpT 1 = k^T / (&irr)R?) (also called the rotational 
diffusion coefficient), with the two stochastic modes de- 
coupled from each other. For a particle propelled with 
velocity V, the direction of motion is itself subject to 
rotational diffusion that leads to a coupling between ro- 
tation and translation. In this case, one can show that 
the (2D projection) mean squared displacement is given 
as: 
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This expression has limiting forms of AL 2 = AD At + 
V 2 At 2 for At < tr and AL 2 = (AD + V 2 TR)At for 
At ^> tr. At short times, the contribution to the dis- 
placement due to the propulsion is linear in time, while 
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FIG. 4: (color online), (a) Velocity determined from the At <C tr behavior for the control (a) and Pt-coated particles 
(o). The solid line is the line of best fit using Eqs. © and ©, with a = 1.2A, A = 5A, fci = 4.4 x 10 11 /xm s , and 
ki = 4.8 x 10 10 fjxaT 2 s~ 1 . (b) Diffusion coefficient determined from the At <C tr behavior for the control and Pt-coated 

particles with the theoretical value ( ) indicated, (c) Rotational diffusion time determined at At S> tr (from the gradient). 

Value at [H2O2]=0 is the theoretical value. 



the Brownian displacement is proportional to the square 
root of time, with these two contributions adding in 
quadrature. At times long compared to the rotational 
diffusion time, rotational diffusion leads to a randomiza- 
tion of the direction of propulsion, and the particle un- 
dergoes a random walk whose step length is the product 
of the propelled velocity V and the rotational diffusion 
time, leading to a substantial enhancement of the effec- 
tive diffusion coefficient over the classical value, namely 
D c ff = D + jV 2 t^. A similar crossover behavior has 
been observed in the motion of tracer beads in the flow 
caused near surfaces covered by bacterial carpets, as well 
as beads propelled by adsorbed bacteria 
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Figure [2] shows the average mean squared displacement 
as a function of elapsed time for half-coated particles im- 
mersed in hydrogen peroxide at various concentrations. 
Each curve in Fig. [2] is produced by averaging over 3000 
frames. In the absence of hydrogen peroxide, the curves 
are linear, indicating the simple diffusive behavior ex- 
pected of a particle undergoing Brownian motion. With 
the addition of hydrogen peroxide they show an increas- 
ingly important parabolic component; the particles are 
being propelled by the local osmotic pressure gradient 
created by the asymmetric chemical reaction. Fitting 
these plots to the limiting forms of the interpolation for- 
mula of Eq. ([I]) yields values for the parameters (V, D, 
and tr) as functions of H2O2 concentration (see Fig. [3]). 
For each concentration, we perform averaging over 20 in- 
dependent movies for the values of the fitting parameters. 

Figure shows the propulsion velocity as a function 
of the hydrogen peroxide concentration. As we expect, as 
the concentration of "fuel" rises the velocity with which 
the Pt half-coated particles are propelled rises from the 
value of zero observed in the absence of fuel; for the con- 
trol particles, without catalyst, no propulsion is observed. 
The propulsion velocity can be calculated using the lat- 
eral gradient of the excess solute concentration C in the 



vicinity of the particle surface. This gradient creates a 
slip velocity V s = (m)hC with the diffusiophoretic mobil- 
ity given as /i = k^TX 2 jr\ where A represents the range 
of the interaction zone between the solute and the par- 
ticle One can then solve the Stokes hydrodynamics 
around the sphere subject to this local slip velocity pat- 
tern on its surface, and find the propulsion velocity. For 
a particle which is half-coated with a material that can 
produce excess particles with diffusion coefficient D a at 
a rate per unit area k, the propulsion velocity is [I2J]: 



V = 



3n u \2 

= — ka\ , 
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(2) 



where the hydrodynamic radius of the solute a is intro- 
duced by way of the Stokes-Einstein relation. 

Since the velocity is directly proportional to the ef- 
fective surface reaction-rate k [Eq. ©], the linear initial 
increase and the subsequent tendency towards saturation 
observed in Fig. 0^ suggests that the Pt-catalyzed break- 
up of H2O2 probably occurs in two stages of formation of 
a Pt(H 2 C>2) complex at a rate per unit area k\ [H20 2 ] vo i, 
in which [H202] V oi has the units of volume percentage, 
followed by a decomposition into water and oxygen at a 
rate per unit area fc 2 . This leads to a Michaelis-Menten 
behavior similar to enzymes [l9j . as opposed to the ki- 
netics of decomposition in homogenous solution, which 
is second order and would produce an upward curvature 
in Fig. [4^. Solving the diffusion-reaction equations with 
the proposed reaction kinetics, we find 
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The solid line in Fig. [4ji is obtained from equations 
© and ©, with a = 1.2 A, A = 5 A, h = 4.4 x 
10 11 #m -2 s -1 , and fc 2 = 4.8 x 10 10 fjm~ 2 s~ l , which are 
reasonable numbers. In contrast to the strong depen- 
dence of the propulsion velocity on hydrogen peroxide, 
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the fitted particle (bare) diffusion coefficient D shown in 
Fig. [4p has little dependence on the presence or absence 
of catalyst and the concentration of hydrogen peroxide, 
and is close to the predicted value from Einstein-Stokes 
relation (a slight residual dependence on hydrogen perox- 
ide concentration may be a result of a coupling between 
these two fitting parameters in the analysis of the dis- 
placement curves). Figure HJ: shows fitted rotational dif- 
fusion time. As mentioned above, the long time effective 
diffusion is substantially enhanced due to the propulsion. 
For the largest hydrogen peroxide concentration, we find 
a value of D c s = 9.0 /xm 2 s _1 , which is an enhancement 
over the purely Brownian value I? by a factor of nearly 
30. 

The rotational diffusion time shows a systematic de- 
crease as a function of hydrogen peroxide concentration 
for the coated particles, with no dependence for the con- 
trols. This could be due to the surface reaction im- 
parting a small net angular velocity uj, and we surmise 
that it might be caused by inadvertent asymmetric Pt 
coverage during the fabrication. Similar to the transla- 
tional diffusion coefficient, we would expect a combina- 
tion Tp 1 + w 2 tr to serve as the renormalized rotational 
diffusion coefficient. This hypothesis is strengthened by 
the observation of cycloid trajectories for some samples at 
high H2O2 concentrations. At first sight this would seem 
an unwelcome side effect — an increase in rotational veloc- 
ity has the effect of decreasing the renormalized diffusion 
coefficient for the long-time behavior of the propelled par- 
ticles. But it also raises the intriguing possibility of de- 
signing a system in which the linear propulsive and rota- 
tional behavior are independently controlled; this would 
constitute a system in which the step size of a random 
walk could be controlled by an external parameter, open- 
ing up the possibility of designing a system capable of 
chemotaxis, in a way analogous to the bacteria E. Coli 
1- 

By contrast to the existing experimental examples of 
autonomous swimmers [H, [3, fiol Il3j |. the propulsion 
mechanism for our platinum/polystyrene particles does 
not involve electrochemical reactions, and thus they real- 
ize a new class of micro- and nano-scale chemical locomo- 
tion. We note that for this class of (self-diffusiophoretic) 
swimmers, spherical geometry is better than a rod geom- 
etry, because the velocity for rods is reduced by a factor 
of the aspect ratio [T3|. Moreover, their true directional 
motion is easier to resolve experimentally as the lower 
friction coefficient of the rods along their length causes 
them to spend more time going along their own direction, 
which in short time probing might be mistaken with di- 
rectional locomotion. 

In conclusion, we show that spatially asymmetric catal- 
ysis at the surface of synthetic micron-scale particles can 
lead to effective autonomous propulsion of spherical col- 
loids. At short times, we observe directed propulsion of 
the particles with velocities in the fim/s range; at longer 



times the direction of motion of the particles is random- 
ized, and the motion of the particles becomes diffusive 
in character with an effective diffusion coefficient which 
is substantially enhanced over the Brownian value. We 
hope that our quantitative characterization of the motion 
of the swimmer sheds some light on the fundamental is- 
sues involved in designing chemical locomotive systems, 
and could inspire new directions for their implementa- 
tion. 
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